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SUMM\RY  PAGE 


THE  PROBLEM 

Failures  in  visual  pursuit,  in  ocular  stabilization  during  head  move¬ 
ment,  and  in  visual  suppression  of  the  vest ibulo-ocular  reflex  (VOR)  are 
characterized  by  saccadic  (high  velocity)  corrections  for  errors  in  the 
smooth  pursuit  velocity  required  to  maintain  foveation  of  targets.  Measures 
of  the  ratio  of  pursuit  eye  velocities  to  velocities  required  by  the  stimu¬ 
lus  situation  together  with  counts  of  saccadic  corrections  are  required  to 
quantitatively  characterize  the  over-all  response  which  is  crucial  to  per¬ 
formance  Ln  many  operational  situations  and  in  neuro-oto logical  examinations 
for  equilibration  disorders. 

FINDINGS 

A  signal  processing  algoritlim  utilizing  available  equipment  was  devel¬ 
oped  to  detect,  count,  and  remove  sac.cades  from  oculomotor  responses  gener¬ 
ated  by  head  oscillation  and  to  provide  estimates  of  the  "response  gain" 
under  several  eondltions  of  visual  stimulation.  After  saccade  removal,  the 
velocity  of  "smooth"  pursuit  or  slow  phase  eye  velocity  may  be  plotted  with 
respect  to  stimulus  velocity  automatically  if  desired.  A  linear  regression 
least  squares  fit  of  the  data  provides  a  quantita'  Ive  estimate  of  the  ' "response 
gain"  (ratio  of  eye  to  stimulus  velocity),  and  the  correlation  coefficient 
provides  an  index  of  the  goodness  of  fit.  This  procedure  may  be  used  to 
assess  the  oculomotor  response  during  either  cyclic  vestibular  or  visual 
s  timulation. 
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INTRODUCTION 


Stabilization  of  visual  fixation  during  natural  head  movement  is 
necessary  for  maintained  high  visual  acuity.  Two  neural  pathways  send 
control  messages  to  the  motor  nuclei  which  drive  the  eyes  smoothly 
during  head  movements;  One  pathway  carrying  information  about  visual 
target  movement  in  the  visual  fields  attempts  to  maintain  foveation  by  a 
"tracking"  or  "smooth  pursuit"  process;  a  second  pathway  originates  in 
the  inner  ear  sensors  of  head  angular  and  linear  acceleration  -  the 
vestibular  organs  -  and  sends  signals  which  are  compensatory  for  the 
resultant  head  movement,  the  vestlbulo-ocular  reflex  (VOR) .  Failure  to 
maintain  foveation  may  produce  a  rapid  compensatory  movement  (saccade) 
of  the  eyes  to  a  new  fixation  point.  These  two  separate  systems  normally 
operate  over  different  ranges  of  target  velocity,  tracking  being  optimum 
from  frequencies  close  to  zero  up  to  1  Hz,  and  the  VOR  being  effective 
from  below  0.1  Hz  to  at  least  4  Hz. 

In  order  to  test  the  function  of  the  VOR,  it  is  common  practice  to 
record  e;'e  movements  electL'lcally  by  detecting  with  skin  surface  electrodes 
the  changing  potentials  produced  across  the  eyes  by  the  moving  corneoretlnal 
potentials  of  the  eyes.  Analytical  determination  of  the  relationship 
between  head  velocity  and  eye  velocity  during  passive  whole-body  rotation 
or  voluntary  head  movements  then  characterizes  the  VOR  (3-5, 7, 9).  TVie 
presence  of  saccades  in  tl\e  ,ye  position  record  complicates  the  analysis 
since  i'  represents  a  different  process  from  that  involved  in  smooth  eye 
movement . 

This  report  describes  a  computer-based  procedure  for  determining, 
in  the  time  domain,  the  output  (eye  velocity)  versus  input  (head  velocity) 
relationship  during  voluntary  head  oscillation,  including  algorithms  to 
detect  and  remove  saccades  before  analysis.  The  procedures,  developed  for 
use  wltli  ciuipmcnt  available  at  tiie  time,  provide  a  basis  for  further  ad¬ 
vances  in  automated  scoring  as  updated  equipment  becomes  available. 

FROCRDURi; 

Head  and  eye  position  signals  were  olitalned  from  htunan  subjects 
according  to  procedures  described  elsewhere  (6),  and  recorded  in  analog 
form  on  nvignetlc  tape,  along  with  calibrations  of  head  and  eye  signal 
channels  for  20  deg  of  left  and  right  movement  with  respect  to  center 
gaze  position. 

Fcir  analysis,  the  magnetic  tape  was  replayed  for  analog-to-dlgital 
( onveralon  and  digital  storage  of  data  on  diskettes.  Analysis  was  then 
carried  out  on  these  diskette  dat.i  flies. 

APPAKA'l'U.S 


Analog  recordings  were  made  on  a  Philips  Mini-Log  4  four-channel 
cassette  recorder.  Recording  mode  was  frequency  modulation  at  a  tape 
speed  of  0.75  Inch  per  sec,  giving  a  signal  channel  bandwidth  of  d.c.  to 
250  Hz.  Lye  and  head  position  signals  were  slmuitaneously  recorded  on 
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paper  chart,  using  a  Hewlett-Packard  Model  775AA  thermal  strip  chart 
recorder. 

The  Hewlett-Packard  Model  3582A  Spectrum  i^lyzer  is  capable  of 
simultaneous  two-channel  analog- to-dlgltal  conversion  and  storage  of 
preselected  time  epochs  of  data  under  remote  control  via  the  standard 
Hewlett-Packard  interface  bus  (HP-IB) .  Tlie  input  circuits  of  the  3582A 
incorporate  appropriate  anti-aliasing  filters  for  the  chosen  sampling 
rate.  Programs  were  written  in  BASIC  to  run  on  a  Hewlett-Packard  9830A 
calculator  which  was  interconnected  via  the  KP-IB  with  the  spectrum 
analyzer.  Under  program  control,  the  spectrum  analyzer  first  was  com¬ 
manded  to  acquire  the  two  channels  of  data  present  at  its  inputs,  gener¬ 
ating  two  corresponding  512-sample  blocks  of  16-blt  digital  data  in  its 
1024  word  memory.  A  typical  time  epoch  was  12.5  sec,  during  which  time 
512  samples  each  of  head  and  eye  position  data  were  obtained,  at  a 
Sampling  rate  of  40.96  samples  per  sec. 

At  the  end  of  each  acquisition  epoch,  the  contents  of  the  spectrum 
analyzer  meiiKiry  were  transferred  via  the  HP-IB  and  under  9830A  control 
to  the  memory  of  tlie  9830A,  and  thence  to  permanent  mass  storage  on  an 
Infotek  Model  FD30A  diskette  unit. 

Analysis  was  done  on  the  y830A  calculator,  using  BASIC  programming. 
D.ita  plots  were  made  by  the  Hewlett-Packard  Model  9862A  X-Y  plotter  and 
printouts  by  the  Hewlett-Packard  Model  987] A  printer,  botli  operating  as 
peripheral  devices  under  the  control  of  the  9830A  calculator. 

MET'HOD 


A  typical  Lime  record  of  head  and  eye  position  data  is  shown  in 
Figure  i  and  will  be  used  as  the  basis  for  this  discussion.  In  this 
case,  voluntary  head  oscillation  at  approximately  1  Hz  over  an  angle  of 
about  40  deg  was  the  forcing  function,  and  the  VOR,  with  partial  suppression 
by  ai tempted  fixation  on  a  spot  of  light  moving  with  the  same  angle  and 
phase  as  the  head,  was  the  response.  Perfect  suppression  would  produce 
a  flat  eye  position  trace,  as  the  eye  would  not  move  in  the  head.  Zero 
suppression  would  result  in  an  eye  trace  matching  tlu^  he. id  waveform,  hut 
180  deg  out  of  phase  with  it.  At  I  Hz,  fixation  was  only  partially  able 
to  suppress  the  VOK  .ind  resulted  in  spor.adic  saccades  to  catch  the 
target . 


'file  objective  was  t.o  rel.ite  eye  velocity  to  head  velocity  throughout 
the  epoi'li,  wltliout  contamination  from  Liie  high  viiloclty  saccades.  The 
s.iccade  detection  and  removal  procedure  are  outlined  in  block  form  in 
Figure  2.  The  operation  was  performed  only  on  the  eye  movement  record 
which  started  out  ;is  the  raw  displacement  dat.i  from  the  diskette  file. 
This  512-point  record  was  diglL.illy  differentiated  by  the  procedure: 


Y 


1  +  L 


.  r 


where  Y  is  the  1th  displacement  sample,  Y  ^  j  is  the  1  +  1st  sample, 
Y'^  is  tlic  ith  velocity  sample,  and  r  is  trie  sampling  rate.  Figure  3 
shows  velocity  data  (head  and  eye)  corresponding  to  Figure  1.  The  eye 
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velocity  record  was  saved  for  future  reraoval  of  saccades,  and  its  duplicate 
was  then  used  for  saccade  detection. 

The  saccade  detection  process  consisted  of  a  digital  high  pass 
filter,  which  removed  ail  but  the  high  velocity  saccades  and  differen¬ 
tiation  noise,  followed  by  full  wave  rectification  and  elimination  of 
subcriterion  peak  velocities.  The  remaining  peaks  were  taken  as  valid 
pointers  to  saccades  in  the  velocity  record.  Implementation  of  the 
digital  high-pass  filter  was  achieved  by  subtracting  from  the  original 
velocity  record  a  low-pass  filtered  version  of  itself.  The  low-pass 
filter  equation  was  the  3-point  moving  average 


^’i  = 


,  +  Y.  +  Y.  ^  , 
1  1  1+1 


where  Y  is  the  ith  sample  of  the  velocity  record  which  is  replaced  in 
the  filtered  record  bv  and  Y  _  j  and  Y,  ^  .  are  the  preceding  and 

next  velocity  samples,  respectively.  Four  iterations  of  this  equation 
were  applied  to  produce  a  low-pass  filter  with  a  cutoff  frequency  of 
approximately  4  Hz,  slope  12  dB/  octave  and  zero  phase  sliift  (8). 

Figure  4  (lower  tracing)  stiows  the  resulting  eye  velocity  record,  which 
was  then  subtracted  from  tlie  record  of  Figure  3  and  full  wave  rectified 
simply  by  the  process  of  taking  absolute  values  of  all  points.  From  the 
full  wave  rectified  record,  as  siiowu  in  Figure  5,  tlie  saccades  could  be 
Located  by  detecting  the  Largest  peak.s  in  eacli  group.  This  process  was 
facilitated  by  first  smootliing  the  record  by  low-jiass  filtering  with  two 
passes  of  the  3-point  moving  average  already  described,  resulting  in  the 
record  shown  in  Figure  6.  The  peaks  of  the  record  were  then  determined 
by  detecting  positive  to  negative  slope  transitions.  All  peaks  below  a 
preassigned  threshold  were  rejected  and  the  indices  of  the  rest  were 
held  In  a  saccade  index  matrix.  At  this  point,  each  saccade  index  was 
used  to  find  the  direction  of  the  corresponding  saccade  by  referring 
back  to  the  record  containing  saccades  before  rectification.  The  sign 
of  the  indexed  sample  in  this  record  Indicated  whether  the  saccade  was 
upgoing  or  dtjwngoing,  and  a  saccade  direction  matrix  was  formed  to 
contain  direction  infonnatlon  for  all  accepted  sccades. 

Sacciide  remov-il  from  the  eye  displacement  or  velocity  racord  was 
performed  by  replacing  five  sample  points  for  each  saccade,  ,is  Indexed 
by  ttie  Saccade  Index  matrix.  A  satisfactory  result  was  obtained  by 
removing  the  indexed  sample,  three  previous  and  one  following;  these 
samples  were  removed  by  replacing  tlieir  values  with  an  unlikely  one,  in 
this  case  tlie  value  999.  Subsequent  processing  always  tested  euch 
sample  for  tlie  value  999,  and  ignored  such  samples.  An  example  of  a 
dispiav’ement  rccon'  with  saccades  removed  is  shown  in  Figure  7,  which 
should  be  compared  with  tlie  rav;  data  of  Figure  1.  Removal  of  saccades 
from  the  velocity  record  is  sliown  in  Figure  8,  wliich  should  be  compared 
w:  til  Figure  3. 

'I'he  true  relationship  Between  eye  velocity  and  head  velocity  was 
tlien  obtained  by  plotting  the  eye  velocity  record  with  saccades  removed 
against  the  corresponding  head  velocity  record.  An  example  of  this  is 


6 


9 


given  in  Figure  9,  where  mov"'  enta  to  the  right  yield  positive  velocltlcb. 
Each  point  represents  a  head  v» loclty  and  the  corresponding  eye  velocity 
at  tliat  point  in  time,  vltli  no  saccades  represented.  The  slope  of  the 
best  fitting  line  for  positive  head  velocities  yields  the  gain  of  the 
VOR  for  head  movements  to  the  right,  and  similarly  tlie  slope  of  the  best 
fitting  line  for  negative  head  velocities  yields  VOR  gain  for  leftward 
lie  ad  movements. 

Best  fitting  lines  by  the  method  of  least  squares  for  botli  leftward 
and  rlglitward  head  movements  liave  been  computed  and  plotted  on  the  data 
of  Figure  9,  and  are  shown  in  Figure  10.  A  line  was  also  computed  and 
plotted  for  tlie  whole  record.  Equations  of  the  best  fitting  lines  and 
regression  coefficients  were  printed  out,  along  with  saccade  number  and 
direction  subtotals. 


result:!  and  DISCU.S.SION 

Comparison  of  Figures  1  anti  7  reveals  t.ie  success  rf  tlie  saccade 
extraction  procedure.  Twenty'-nlne  s.iccades  were  detected  and  removed 
from  the  data  of  Figaro  1,  le.ivlng  a  displacement  recortl  consisting 
prlnvu  lly  of  the  slow  compi'i.ont  of  the  VOR.  When  the  same  saccades  were 
removed  from  the  velocity  tlata  of  Figure  3,  the  result  shown  In  Figure  8 
was  ohtalntul.  Eye  velocity  was  then  plotted  /igaliisL  head  velocity,  and 
the  result  Is  shown  lu  Figure  9,  where  each  ilot  represents  one  eye 
velocity  sample  and  the  head  velocity  tairrespondl ng  to  It.  Most  of  the 
dots  are  clustered  about  lines  through  the  origin  In  the  second  and 
fourth  ((uadrants,  and  would  all  fall  on  the  bisectors  If  the  gain  was 
unit)  with  zero  phase  shift.  The  actual  gain  was  determiited  by  fitting 
a  regression  line,  l)v  'he  metliod  of  least  squares,  to  tlie  data  In  each 
side  of  the  vertical  axis;  dat.i  on  left  side  represent^  head  movements 
Low.'trd  the  left  while  those  on  the  right  side  are  representative  of  head 
movements  to  the  right.  Results  of  the  Linear  regression  for  these  data 
are  shown  In  Figure  10.  A  slightly  higher  gain  with  head  movements  to 
Llie  right  is  Indicated  (0.7/7  vs.  0.671),  although  tlie  spread  is  higher 
with  leftward  head  inovemeuts,  A  regression  line  fitted  to  all  data 
jiroduces  an  1  ntermed  I /i  t  e  value  ol  g.ilii  (0.689)  with  better  fit  (t  "  0.886) 
than  either  d  1  rei' '•  loii.i I  line. 

lu  order  to  assess  t  Ik'  value  of  s;u-c;ide  extr.ict  loii  In  deterniln  1  iig 
these  p.'irameLers,  linear  regi  ess  Ions  were  computed  for  the  eye  veloivlty 
record  without  removal  of  sici.ides.  These  data  are  shown  In  Figure  11, 
where  the  dots  h.ivc  lieen  connect  c'd  with  llin's  to  show  t  lu  saccade  tra)ect,i)ry 
more  clearly.  Regression  analysis  produeed  unns.ihle  }’,aln  data  In  this 
c.ise,  with  slopes  skewml  by  high  saec.ide  velocities,  .ind  low  coefficients 
of  correlation  as  a  result  of  the  spre.id. 

The  procedures  deserlhed  have  thufi  been  shown  to  be  success  fnl  lu 
prodiu'lng  .i  meaningful  analysis  of  eye  versus  head  velocity  r(*  I  ation.shl  ps 
under  tlie  circunisLanees  given.  Two  aspects  of  the  procedures  are  empirical 
and  could  give  rise  to  less  optimum  analyses  under  certain  other  cl rcumstances 
These  are  the  value  of  tlie  threshold  level  for  peaks  used  to  identify 
valid  saccades  (c'  shed  line  in  Figure  6),  and  the  number  of  samples 
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removed  before  and  after  tlve  peak  index  from  the  velocity  record.  The 
fi'rmer  1b  related  to  .saccade  velocity  which,  under  unchanging  conditions 
of  Input  to  the  VOR  (head  movement),  might  be  expected  to  be  ’-easonably 
const. iOt  for  each  saccade  as  long  as  it  occur.s  at  the  same  he. id  velocity 
in  each  cycle,  ^ny  head  velocity  or  eye  amplitude  dependence  (2)  of  the 
saccade  velocity  may  cause  error.';  to  occur  In  which  true  saccades  are 
not  eliminated  because  their  velocities  are  to  ■  low  for  ttie  chosen 
threshold.  The  effectiveness  of  saccade.  removal  In  producing  a  close- 
fitting  regression  line  is  also  related  to  how  well  each  saccade  is 
ts'inoved  once  It  is  detected  by  the  program,  and  this  bolls  down  to  how 
many  sampies  before  and  liow  many  samples  after  the  peak  Index  are  eliminated. 
Clearly,  If  some  samples  during  the  saccade  are  left  In,  a  few  straggling 
points  wllJ  remain  to  lessen  the  regression  line  goodness  of  fit.  Such 
Is  obviously  tlie  case  in  Figure  10.  Alternatively,  if  too  m.any  samples 
.are  eliminated  before  and/or  after  each  saccade  ii\dex,  valid  points  will 
be  lost  from  tl\e  .analysis.  This  could  be  of  Import.ance  in  evaluating 
the  Vt)R  in  vcstibui.ar  disorders  which  may  Involve  high  velocities  in 
particular  (1)  since  the  saccades  predominantly  occur  during  the  highest 
velocities  of  head  nuivement.  Careful  assessment  of  this  factor  is  therefore 
necessary  tn  the  application  of  the  technique  reported  here. 
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